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^ I Abstract 

■ A comparative study of various parametrizations of the dark energy equation of 
^ I state is made. Astrophysical constraints from LSS, CMB and BBN are laid down to 

test the physical viability and cosmological compatibility of these parametrizations. 
. A critical evaluation of the 4-index parametrizations reveals that Hannestad-Mortsell 

i as well as Lee parametrizations are simple and transparent in probing the evolution 

' of the dark energy during the expansion history of the universe and they satisfy the 

■ LSS, CMB and BBN constraints on the dark energy density parameter for the best fit 
r~| . values. 

S ; 1 Introduction 

■ 

^ . The recently discovered 16 Type la supernovae (SN la) with the Hubble Telescope by Riess 
• 1— I . et al. P provide a distinct scenario and conclusive evidence of the decelerating universe in 
^ i the past (2; > 1) evolving into the present day accelerating universe. Thus, the existence 
I of dark energy, which accelerates the cosmic expansion, has been firmly established and the 
present magnitude of its energy density has been precisely determined 0. The focus is 
now on the evolution of dark energy and its equation of state with the cosmic expansion. 
The simplest and the most natural candidate for dark energy is the cosmological constant 
A with a constant energy density p\ and a fixed equation of state parameter w = —1. But 
the recent analysis of the SN la data [3 E] indicates that the time varying dark energy 
gives a better fit to observational data than a cosmological constant which has only one 
free parameter fi^j. Assuming the time varying equation of state, an important issue arises 
whether the dynamical dark energy can cross the phantom divide line (PDL) w = —1 during 
evolution? To put it more explicitly, can quintessence transform into phantom dark energy 
during the course of its evolution? There are divergent views on this question [5-8] but it is 
found that such parametrizations which cross the PDL seem to provide slightly better fit to 
observational data [H]. 
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As regards the rapidity of the evolution of dark energy, there are two contrasting views at 
present. Riess et ah pOj and Jassal et al. have argued that the current SN la observational 
data is inconsistent with the rapid evolution of w{z). On the other hand, Bassett et al. ^Hj 
and Corasaniti et al. ^I] contend that the inconsistency arises on account of the inadequacy 
of 2-index parametrization.. They claim that the rapid variation in w{z) in fact provides 
a better fit to the SN la 'gold set' observations, even after including CMB and large scale 
structure data [llj. According to our analysis, the rapidity of variation |^|, apart from 
other factors, depends on the absolute value of \w\. For quintessence models —l<w<—^, 
whereas w has no lower bound for phantom models. Hence, the equation of state of dark 
energy varies more rapidly for phantom models with large \w\ than for quintessence models, 
other factors remaining unchanged. The fact that the observational data by Caldwell et al. 
[T^ including SN la and galaxy clustering, shows a bias towards phantom models might 
explain why rapid variation in w{z) provides a better fit to observational data. 

It is difficult to know the exact functional form of the equation of state parameter w{z) 
observationally as such different parametric forms namely one index parametrization by Gong 
et al. jni, two index parametrization [9, 14-21], four index parametrization by Hannestad 
et al. [22] and Lee ^H] , power law parametrization ^3 |211 123 ) more complicated forms 
[26-29] and non-parametric forms Bassett et al. ^U], Corasaniti et al. fUEn] and Sahni et 
al. jSI] of w{z) have been used to simulate the behavior of evolving equation of state of dark 
energy and comparison made with the SN la observations. Alternatively, the dark energy 
density also can be parametrized as power law expansion [31-36] or piece wise constant 
parametrization j21 IH EH IHH] • The SN la observations essentially measure the luminosity 
distance diiz) which when compared with the theoretical parametric values, yields the best 
fit values of the parameters. But it may be stressed that these observational tests alone 
are not adequate to distinguish between the cosmological constant and the dynamical dark 
energy. If we start with the assumption w = constant, the observational data favors w = —1. 
In case of dynamical dark energy, the current observational data does not yield any conclusive 
constraint over the equation of state parameter w{z). In fact, large downward deviations of 
the equation of state parameter from w = —1 favor phantom dark energy. Therefore, we need 
to devise some new tests for probing dark energy. Odintsov and Nojiri have discussed 
the thermodynamics of the evolution of dark energy and have stressed that the phantom 
stage might be a transient period in cosmic evolution. So far we have been exploring the 
properties of dark energy (its equation of state and the density parameter fix) using only 
cosmological observations. A new approach might be to use astrophysical constraints to 
probe the various dark energy parametrizations of w{z) as discussed in this paper. 

For making a comparative study of various parametrizations and their viability, we need 
some criterion. In our previous papers [40-42], we introduced the concept of 'integrated 
tracking' and outlined certain astrophysical constraints to be satisfied by the dark energy 
fields for realistic tracking. Recently, Upadhye et al. 20j have discussed the dark energy 
constraint from CMB data at decoupling epoch z = 1100. The astrophysical constraints 
discussed below arise from well established cosmological observations at different redshifts. 

(A) LSS Constraint 

During the galaxy formation era (1 < z < 3), dark energy density must be sub- 
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dominant to matter density; accordingly Qx < 0.5 [IS] . 

(B) CMB Constraint 

Caldwell ^^ and Upadhye et al. jSH] have argued that {ilx)dec < 0.1 at z = 1100. It 
imposes a high redshift upper bound fix for any viable parametrization of the equation of 
state parameter w{z). 

(C) BEN Constraint 

The presence of dark energy until nucleosynthesis epoch, should not disturb the observed 
Helium abundance in the universe which is regarded as the foundation stone of the Big Bang 
Theory. According to Johri |4rj {flx)BBN < 0.14 at 2; = 10^° whereas the latest analysis 
of Cyburt jlH] constrains {flx)BBN < 0.21 at z = 10^°. This puts a stringent high redshift 
limit on dark energy. 

In this paper, we have discussed the cosmological compatibility and the range of validity of 
six dark energy parametrizations, the first four involving 2-parameters and the fifth and sixth 
involving 4-parameters. More parameters mean more degrees of freedom for adaptability to 
observations, at the same time more degeneracies in the determination of parameters. All 
the parametrizations are critically examined on the basis of the above constraints to test how 
faithfully they represent the behavior of dynamical dark energy in the expanding universe. 

The outline of the paper is as follows. In section |2l we have discussed the expansion 
history of the universe. Assuming a spatially fiat universe, the field equations give the Hub- 
ble expansion and deceleration parameter in terms of function f{z) which involves integral 
of the varying equation of state parameter of dark energy. In section 3, we have discussed 
2- index parametrization models namely linear red shift parametrization jEl ^] , Chevallier- 
Polarski-Linder parametrization fW\ [T7] . Jassal-Bagla-Padmanabhan Parametrization jH] 
and Upadhye-Ishak-Steinhardt parametrization ^20j. In section 4, we have critically ex- 
amined Hannestad-Mortsell parametrization [22] and Lee's model j2S| which involve 4 pa- 
rameters imparting four degrees of freedom to choose them. In section 5, we use interpolation 
technique to study the behavior of the dark energy under the assumption of slowly varying 
equation of state. In section 6, we conclude with some remarks on parametrization methods. 



2 Expansion history of the universe 

Assuming a spatially fiat (k = 0) Friedmann universe, the field equations are 

H' = H'^[ni,{i + zf + ny{z)] (i) 

and 

= w{z)nx (2) 

where if = ^ is the Hubble constant, q = is the deceleration parameter, 'w{z) = 

is the equation of state of dark energy and VLx = 1 — ^Im is the cosmic dark energy density 
parameter. The dark energy density parameter is defined as 

flx = ^— (3) 

PC 
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with 



and 



Px 



exp 



l+wiz') 



dz 



(4) 



(5) 



10 1 + ^' 

We can test any parametrization by taking the best fit values to the observed luminosity 
distance of SN la given by 



and applying constraints (A), (B) and (C) to check its compatibility with the cosmological 
observations. In Table 1, we have tabulated the dark energy parameters Vt^ at LSS formation, 
decoupling and BBN corresponding to ^ = 1 — 3, 1100 and lO^'' respectively. The transition 
redshift may also be derived from the best fit values. 



3 Two index parametrization 

I. Linear-redshift parametrization 

The equation of state parameter w{z) is given by Huterer et al. 
dS] as 

w{z) = Wo + w'z, 
Inserting Eq. (7) for w{z) in Eq.(3), we get 

i ' A/I 



w 



/dw\ 
\dz J z=o 



and Weller et al. 



(7) 



1 + 



l + z) 



-3{wo—w') 



exp{—3w'z) 



It has been used by Riess et al. [T] for probing SN la observations at 2; < 1. The best fit 
values to SN la 'gold set' data [iHj are Wo = —1.4, w' = 1.67 and fl%j = 0.30. Hence, the 
parameters favor dark energy of phantom origin. In the distant past { z » 1) , w ^ 00 and 
fix 1 and in the distant future {1 + z ^ 0) , w = Wq — w' = —3.07 and Qx — ^ 1. At the 
transition red shift zt = 0.39, we get {flx)T = 0.443. The parametrization is not viable as 
it diverges for z » 1 and therefore incompatible with the constraints (B) and (C). 

II. Chevallier-Polarski-Linder parametrization 

The equation of state parameter w{z) is given as [121 HZj 

W{z) = Wo + Wi— — (9) 

l + z 

It indicates rapid variation of w{z) at z = which goes on 



dw 



On differentiation, — jj+zf 
decreasing with increasing z. Eq.(3) gives 



n 



^^1 + z' 



-I -1 



(10) 
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The best fit values of Wo-, wi and to SN la 'gold set' data [29l Eoj turn out to be —1.6, 
3.3 and 0.3 respectively. At z = 0, one has w{z) = wq = —1.6. and the parameters suggest 
that the dark energy is of phantom origin. For z » 1, w ^ wq + Wi = 1.7 and fix 1- 
In the distant future {1 + z ^ 0), w ^ —oo and fix — ^ 1- The dark energy parameter 
{flx)T = 0.452 at the transition red shift zt = 0.35. The CMB and BBN constraints for fix 
are violated. Hence, it is incompatible with observation. 

III. Jassal-Bagla-Padmanabhan parametrization 

The equation of state parameter w{z) is given by Jassal et al. 



as 



As ^ 

az 

gives 



■w\(l-z) 



wiz. 



w(z) 



increases from z 



Wo + 



WiZ 



to z = 1 thereafter starts decreasing. Eq.(3) 



fl 



x[z) 



1 + 



fp 



l + z) 



-3 Wo 



exp 



Swiz"^ 



2(1+^)2 

-1.9, wi = 



(12) 



The best fit values to SN la 'gold set' data are Wo = —l.9,wi = 6.6 and fl%.j = 0.3. One 
has w{z) = wq = —1.9 at z = 0. Thus, the parameters suggest that the dark energy is of 
phantom origin. In the distant past for z » 1, w{z) —>■ wq = —1.9 and fix 0. In the 
distant future {1 + z —>■ 0) , w —>■ —oc and fix ^ 1- At zt = 0.3, one gets {flx)T = 0.467. 
Further, it satisfies all the three observational constraints (A), (B) and (C). 

IV. Upadhye-Ishak-Steinhardt parametrization 

Upadhye et al. have parametrized w{z) as |20j 



w{z) 



Wo + WiZ 
Wo + Wi 



z < 1 

z > 1 



(13) 



The dark energy density contrast is given by 



flxiz) 



1 + ^(1 + ^)-3{«'o-«'i) exp(-3wi^) 



1 + ^(1 + ^)-3(«'0+wi)g-3ti;i(l-21n2) 



if z<\ 
%i z^l 



(14) 



The best fit values of parameters for SN la 'gold set', galaxy power spectrum and CMB 
power spectrum data are wq = —1.38, wi = 1.2 and fl^^^ = 0.31. The transition red-shift zt 
turns out to be 0.44. At z = 0, one has w = wq = —1.38 and this suggests that the dark 
energy has phantom origin. In the distant past [z » 1) , w —0.18 and fix ~* and in the 
distant future (1 + z — > 0), w = i^o — u^i = —2.58 and fix 1- One gets {flx)T = 0.392 at 
transition redshift. The parameterization also satisfies the LSS, CMB and BBN constraints. 



4 Four index parametrization 

I. Hannestad-Mortsell parametrization 
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Let us now consider Hannestad parametrization which involves 4-parameters. 



w{z) 



WoWi 



WiOP + WoO^ 



■ 1 + Zs 



w- 



+ w- 



-1/ 1+z \P 



(15) 



where Wq and Wi are the asymptotic values of w{z) in the distant future (1 + z — >• 0) and in 
the distant past {z — > oo) respectively. The a<j and p are the scale factor at the change over 
and the duration of the change over in w respectively. 

It follows that the equation of state at the present epoch {z = 0) is given by 



w 



Wi + WoOFs 

Differentiating Eq.(15) with respect to z, we get 



(16) 



w 



dz 



{wi - Wo)p{l + zy-^al 
[1 + {as/ a)P] [wi + Wo{as/a) 



(17) 



Hence, the equation of state parameter of dark energy varies with the redshift unless Wi = Wo 
or w{z) takes asymptotic values l + 2;-^ooorl + z— i^O. The variation gradient |^| varies 
directly as \w\, \wi — Wo\, p and inversely as af(l + zy. 
Inserting for w{z) from Eq.(15) into Eq.(3), we get 



n 



x[z) 



-3uii 



wi + WoaP {1 + zf 

{Wi + WoO^s) (1 + z) 



-3(u)o_tui) ' 



n -1 



Taking the best fit values for the combined CMB, LSS and SN la data [22j, Wo = —1.8, wi = 
—0.4 , p = 3.41 and = 0.50 with a prior = 0.38, one has w{z = 0) = —1.38 at z = 0. 
So, the parameters suggest that the dark energy is of phantom origin. In the distant past ( 
z » 1), w{z) —>■ wi = —0.4 and Qx 0. In the distant future {1 + z —>■ 0) , w ^ wq = —1.8 
and Qx ~^ 1- At transition redshift zt = 0.39, the dark energy parameter {Qx)t = 0.333. 
Further, all the observational constraints are satisfied by this parametrization. 

II. Lee parametrization 

The equation of state parameter w{z) is parametrized by Lee as 



w{z) = w. 



Wo exp{px) + exp{pxc 
exp(px) + exp{pxc) 



where 



a: = In a = — ln(l + z) 



(19) 
(20) 
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and the symbols Wr and Xc have the usual meaning as in 
fix(^) is given by 



The dark energy parameter 



1 



0° 



a + a 



eq) 



1 + 



-4/p" 



The parameters is chosen to be 1/3 using the tracking condition and is taken as —3. 
The other parameters are obtained by analyzing the separation of CMB peaks and the time 
variation of the fine structure constant. For = 0.27 and Xc = —2.64, the range of p is 
taken as 1.5 < j9 < 3.9. At z = 0, one has —0.975 < w{z = 0) < —1.0 corresponding to 
1.5 < p < 3.9. Hence, the parameters suggest that the dark energy is of quintessence origin. 
For p = 1.5 and 3.9, the dark energy parameter Qx is 0.353 and 0.333 at the transition red 
shift zt = 0.74 — 0.76 respectively. In the distant past for z » 1, w{z) Wr = 0.333 and 



1 + 



In the distant future, 



w 



WrWQ 



■1 and 



1. In this case, all 



the three observational constraints (A), (B) and (C) are satisfied. 



5 Interpolation of w(z) for slowly varying equation of 
state 

If we go by the analysis of Riess et al. PP and Jassal et al. [5], the current SN la data is 
inconsistent with rapid evolution of dark energy. Therefore, we can use 'integrated tracking' 
and interpolation techniques applicable to slow time varying equation of state |41j. It was 
shown |12] that the scalar fields with slowly varying equation of state which satisfy tracking 
criteria, are compatible with astrophysical constraints (A) and (C) outlined under 'integrated 
tracking'. According to our analysis jH], based on integrated tracking, we have 

nx < 0.14, w ~ -0.035 at z = 10^° 
nx = 0.66, w = -0.77 at z = 
nx = 0.5, w = -0.66 at zt = 0.414 

The above data is consistent with the transition redshift given by Riess et al. llj and 
{nx)BBN given by Cyburt et al. 



6 Conclusions 

We have investigated the behavior of dynamical dark energy by taking parametric represen- 
tations of the equation of state. Out of the 2-parameter models of the equation of state, 
Jassal-Bagala-Padmanabhan parametrization and Upadhye-Ishak-Steinhardt parametriza- 
tions are found to be valid and compatible with the astrophysical constraints over a wide 
range of redshift. 

Bassett et al. and Corasaniti et al. [HIl E] have discussed the limitations of 2-index 
parametrizations. Most of them track well at low redshifts z < 0.2 but if we explore dark 
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energy for regions z > 1 and w << —1 (beyond quintessence models), very rapid variation in 
w{z) can be found. Higher order parametrizations are more suitable for probing the nature of 
dark energy and its evolution since more parameters give more freedom to fit in observational 
data but at the same time it gives rise to more degeneracies in the determination of the 
parameters, 

The Hannestad-Mortsell model and Lee 4-parameter model of the equation of state pro- 
vide well- behaved representations of dark energy evolution over a long range of redshift; 
they admit rapid variation of the equation of state as well. Since, there is no lower bound 
on w{z) for the phantom models, ^lx{z) also decreases steeply in the early universe if 1 + wi 
is negative. 
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Table 1: The dark energy parameter fix at LSS formation, decoupling and BBN corre- 
ponding to 2; = 1 — 3, 1100 and 10^° respectively. 



Models 


Reference 




Q 








z = l-3 


1100 




Linear red shift 




0.371-0.957 


1.0 


1.0 


Chevallier-Polarski-Linder 


[mini 


0.362-0.621 


1.0 


1.0 


Jassal-Bagala-Padmanabhan 


M 


0.348-0.184 


0.0 


0.0 


Upadhye-Ishak-Steinhardt 


1201 


0.276-0.208 


0.012 


0.0 


Hannestad-Mortsell 


1221 


0.170-0.066 


0.0 


0.0 


Lee (p = 1.5) 


ESI 


0.270-0.055 


0.052 


0.177 


Lee (p = 3.9) 


m 


0.253-0.041 


0.054 


0.185 
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